NASA TECHNICAL NOTE 






f63f(c 

NASA TN D-6927 


uo 

■< 



C O y 




MEASURED AND CALCULATED MEAN-FLOW 
PROPERTIES OF A TWO-DIMENSIONAL, 
HYPERSONIC, TURBULENT WAKE 


by Richard D. Wagner 

Langley Research Center 
Hampton, Va. 23365 

I 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION . WASHINGTON, D. C. • NOVEMBER 1972 



1. Report No. 2. Government Accession No. 

NASA TN D-6927 

3. Recipient's Catalog No. 

4. Title and Subtitle 

MEASURED AND CALCULATED MEAN- FLOW PROPERTIES 
OF A TWO-DIMENSIONAL, HYPERSONIC, TURBULENT WAKE 

5. Report Date 

November 1972 

6. Performing Organization Code 

7. Author(s) 

Richard D. Wagner 

8. Performing Organization Report No. 

L-8158 

10. Work Unit No. 

501-06-08-01 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, Va. 23365 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

The hypersonic turbulent wake produced by a wedge was studied experimentally and its 
properties were compared with predictions obtained from a numerical computation procedure. 
In the computation procedure several models for the eddy-viscosity formulation of the turbu- 
lent transport were examined. Conventional-defect models and a modified mixing-length 
model were found to yield good predictions of the experimental data. The classical mixing- 
length model gave unrealistic results. The experimental data displayed similarity when 
velocity and temperature defects were scaled by the maximum defects and the transverse 
coordinate was scaled by the velocity-defect half-width. 


17. Key Words (Suggested by Author(s)) 

Hypersonic turbulent wake 
Experimental data 
Numerical computations 
Shear flow 
Turbulent flow 


18. Distribution Statement 
Unclassified ■ 

- Unlimited 


19. Security Oassif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified 

Unclassified 

33 

$3.00 


For sale by the National Technical Information Service, Springfield, Virginia 22151 






















MEASURED AND CALCULATED MEAN- FLOW PROPERTIES OF A 
TWO-DIMENSIONAL, HYPERSONIC, TURBULENT WAKE 

By Richard D. Wagner 
Langley Research Center 

SUMMARY 

The hypersonic turbulent wake produced by a wedge was studied experimentally and 
its properties were compared with predictions obtained from a numerical computation pro- 
cedure. In the computation procedure several models for the eddy -viscosity formulation 
of the turbulent transport were examined. Conventional-defect models and a modified 
mixing-length model were found to yield good predictions of the experimental data. The 
classical mixing-length model gave unrealistic results. The experimental data displayed 
similarity when velocity and temperature defects were scaled by the maximum defects 
and the transverse coordinate was scaled by the velocity -defect half -width. 

INTRODUCTION 

Current interest in hypersonic turbulent shear flows stems from several practical 
problems requiring knowledge of the flow properties in such free shear layers. One such 
problem is the use of vehicle wake observables for entry-vehicle discrimination (ref. 1); 
another is determining the effects of impingement of a free shear layer on the surface of 
a hypersonic vehicle. These shear layers are often produced by the interaction of shocks 
from components of the vehicle (ref. 2). Each of these problems is characterized by a 
close coupling of the developing time-averaged mean flow and the flow turbulence. 

At present, numerical techniques exist which compute turbulent shear flows with 
phenomenological models for the turbulent transport. Usually, these models for hyper- 
sonic turbulent free shear layers are seemingly logical extensions of models which have 
given reasonable results in comparisons with low-speed experimental data. However, 
experiments are needed to assess the applicability of these extensions to the hypersonic 
turbulent shear layer. 

The present paper reports the results of experiments conducted in the Langley 
22 -inch helium tunnel to determine flow characteristics in the turbulent wake produced 
by a 45° total-angle wedge. The measurements made include tests to ascertain the tur- 
bulent character of the wake, and surveys of the wake pitot pressure, total temperature. 



and static pressure. Experimental data at an initial station are used to start computa- 
tions of the subsequent wake development by inserting various turbulent transport models 
into a numerical computation procedure. The data downstream of the initial station are 
compared with computed distributions of flow properties at several stations in the wake. 

SYMBOLS 


C = PM 

Cy 

d 

h 

i 

k 

I 

h 

M 

P 


Por 

Pt 

Pt,2 

R 


constant in intermittency distribution 
specific heat at constant pressure 
diameter of static probe 
static enthalpy 

equal to 0 for two-dimensional flow and 1 for axisymmetric flow 
constants in eddy -viscosity formulation (eqs. (12) and (13)) 
thermal conductivity 
mixing length 

constant in mixing-length formulation 

Mach number 

static pressure 

orifice static pressure 

stagnation pressure 

pitot pressure 

local Reynolds number , 


2 



Rt 

Rtr 


T 



t 

u,v 

x,y,z 

yi/2 

y 

yx 

5 

6 * 

Ij 

e 

M 

p 


turbulent Reynolds number 
transition Reynolds number 

put 

free-stream Reynolds number based on wedge thickness, °° °° 

P oo 

characteristic body dimension; rjj = t for two-dimensional body 

static temperature 

adiabatic wall temperature 

total temperature 

wedge base thickness 

velocity components in x- and y-directions 

Cartesian coordinates (fig. 1) 

half-width for velocity defect 

ratio of specific heats 

intermittency factor 

wake width 

displacement thickness 

transformed coordinates (eqs. (6)) 

momentum thickness 

wedge half-angle 

viscosity 

density 



CT Prandtl number 

T Von Mises variable 

Subscripts: 

e local wake -edge value 

i initial station 

I local value in wake, or laminar quantity 

o value at y = 0 

sd sticking distance 

T turbulent quantity 


tr 

at transition location 

x,y 

X 

or 

y partial derivative 

x,y 

X 

or 

y partial derivative 


1 

or 

T partial derivative 

00 

free- 

stream value 


A bar over a symbol denotes a nondimensional quantity. 

APPARATUS AND TESTS 
Test Facility 

The experiments were conducted in the Langley 22 -inch helium tunnel. This 
blowdown-type tunnel has a contoured nozzle and a 0.559-meter-diameter test section. 
The free -stream Mach number in the test core is nearly constant at any given stagnation 
pressure (i.e., the axial and lateral Mach number gradients are negligible). A calibra- 
tion and description of the facility are given in reference 3. For the present tests the 
stagnation pressure and temperature were 6.89 x 106 N/m2 and 370 K, respectively. The 
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corresponding test-core Mach number was 20.4 and the unit Reynolds number was 
1.34 X per meter. 


Model and Instrumentation 

The wake was generated by a 45° total-angle wedge which had a base thickness of 
2.54 cm and spanned the tunnel test section, (See fig. 1.) The wedge was constructed of 
fiber glass (Hetron 72 laminate) with an inlaid, stainless-steel, sharp leading edge. Five 
static -pres sure orifices (0.23 cm in diameter) were placed in the base to measure base 
pressure. 

A splitter plate could be attached to the base of the wedge to form a rearward- 
facing step with a depth equal to one-half of the base thickness. The plate was 0.46 meter 
long and spanned the base of the wedge. Fifteen static -pressure orifices, 0.23 cm in 
diameter, were placed along the midspan line of the plate and connected to capacitance - 
type pressure transducers. 

Before each test the wedge surface was at room temperature. The stagnation tem- 
perature, 370 K, was selected so that room temperature would be near the adiabatic wall 
temperature for the flow over the windward surface of the wedge. A recovery factor at 
the model surface equal to the square root of the Prandtl number was assumed since the 
flow on the model was laminar; for the wedge-flow Mach number, 2.84, this value is near 
the exact value given by laminar theory, (See ref. 4.) Theoretical analyses for the 
recovery factor of the base flow are not available, and the adiabatic wall temperature on 
the base of the wedge may not equal room temperature. Nevertheless, the low heat- 
transfer coefficients with the small temperature potentials that are likely to exist should 
lead to negligible heat transfer in the base region. Thus, the present flow conditions may 
be assumed to be essentially adiabatic. 

Surveys of wake flow properties were made at several chordwise stations behind 
the wedge. The survey rake is sketched in figure 1. The pitot tube was a 0.152-cm- 
diameter tube which was flattened to 0.102 cm across the flats; the tube was connected to 
a diaphragm-type pressure transducer located on the rake support inside the tunnel test 
section. The total-temperature probe was a shielded chromel-alumel thermocouple with 
an outside diameter of 0.094 cm. The third probe shown in figure 1 was a hot-wire probe 
connected to a constant- current anemometer with a frequency response of 100 kHz. The 
wire was tungsten, 3.56 micrometers in diameter and about 0.1 cm long. The hot wire 
was operated at a fixed current during the surveys, and only the rms of the voltage fluc- 
tuations from the compensating amplifier was monitored. The hot wire was used only to 
assist in detecting the character of the wake (i.e., whether laminar or turbulent). 
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During the 50-second run time available, the survey rake was moved across the 
wake at about 0.17 cm/sec. The probes started about 1 cm below the wake plane of sym- 
metry (confirmed by the pitot-pressure symmetry) and traversed to about 7.6 cm above 
the plane of symmetry. The pitot tube moved in a plane passing through the model mid- 
span, and the total-temperature and hot-wire probes translated in parallel planes, dis- 
placed 0.95 cm on either side of the pitot tube. Tests at fixed positions indicated that 
the speed of the probe response was sufficient to prevent errors due to the continuous 
probe movement. 

Measurements of the low static pressures in the wake required the full test time 
available; hence, detailed surveys were not practical, and only the static pressure along 
the wake plane of symmetry was measured. The static-pressure probe was a 0.317-cm- 
diameter tube with a 17.7° half-angle cone tip. (See fig. 2.) A static -pressure orifice 
was located 5.33 cm from the cone tip. Inviscid pressure ratios (obtained from calcula- 
tions by the inviscid method of characteristics) at the orifice as a function of Mach num- 
ber are shown also in figure 2. For the present tests the local Mach number was always 
less than about 6, so that the inviscid pressure would be slightly less than the local static 
pressure ahead of the probe. The viscous effects on a static-pressure probe having a 
geometry similar to the present probe have been studied by Behrens. (See ref. 5.) 
Although Behrens' experiments on a 10*^ half-angle cone followed by a cylinder were in 
air, neither the effect of the present different forebody (17.7° half-angle cone) nor the 
effect of the different specific-heat ratios is believed to be large; thus, the calibration 
reported by Behrens should give a good estimate of the viscous effects on the present 
static-pressure probe. This estimate implies that the maximum viscous correction 
would be an increase of about 5 percent. That is, the measured pressure would be 5 per- 
cent greater than the inviscid pressure. In view of these offsetting differences, the pres- 
ent uncorrected static pressure might be expected to be close to the true static pressure 
ahead of the probe; further data to support this contention will be discussed subsequently. 

RESULTS AND DISCUSSION 
Survey Data 

Schlieren and electron-beam photographs of the wake are shown in figure 3, The 
flow features, as identified in figure 3(c), are typical of those reported elsewhere in the 
literature. (See ref. 6.) Downstream of the base recirculation region and the wake neck, 
the flow is divided into three identifiable regions; The central viscous wake, the inviscid 
flow surrounding the viscous wake, and the inviscid flow above the wake shock. 

Pitot-pressure and total-temperature profiles are shown in figures 4 and 5 for 
10 positions behind the wedge. (The total-temperature probe was inoperative for the 
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survey at x/t = 4.) The profiles extended through the upper half of the viscous wake, 
the wake recompression shock, and part of the inviscid flow above the wake recompres- 
sion shock. The spreading of the viscous wake and wake recompression shock is evident 
in both the pitot-pressure and total-temperature profiles. At the recompression shock, 
the total-temperature profiles show a decrease in T^ which could be due to a transient 
caused by the discontinuity in recovery temperature of the probe shield as it crosses the 
shock. Everywhere, the total temperature is either less than or equal to the free-stream 
total temperature. This effect is not understood since some increase above free-stream 
total temperature would be expected in view of the overshoot that would occur in the wedge 
boundary layer for adiabatic conditions. 

The surveys did not extend all the way to the wedge bow shock since this flow region 
was not of primary interest; this flow should be simply that resulting from the centered 
expansion of the wedge flow and the interaction with the expansion waves reflected from 
the bow wave. The centered expansion of the wedge flow was verified by determining 
lines of constant pitot pressure in the flow above the recompression shock. (See fig. 6.) 
The center of the expansion (obtained by averaging over the intersections of pairs of lin- 
ear curves fitted to the data by least-squares error) is located about four boundary-layer 
displacement thicknesses upstream of the base and slightly above the wedge boundary- 
layer edge. Theoretical analysis (ref. 7) gives a comparable distance for the upstream 
effect of a corner on an approaching laminar boundary layer. The upstream boundary 
layer "feels” the corner by propagation through the subsonic portion of the boundary layer. 
It should be noted that some error is introduced in the location of the expansion center 
since the region shown in figure 6 includes portions of the expansion fan influenced by 
bow-wave reflections. The limiting characteristic, the first wave of the expansion fan 
that is reflected from the bow wave, is shown to cross the lines of constant pitot pressure. 
(The limiting characteristic and the bow shock were calculated by an inviscid method of 
characteristics.) However, the local Mach number varies from about 6 on the first line 
of constant pitot pressure shown to about 10 on the last. At these high Mach numbers, 
one would not expect significant bending of the expansion-fan rays by the reflections. 

The static-pressure distribution obtained along the wake plane of symmetry with 
the static-pressure probe is shown in figure 7 along with static pressures measured on a 
splitter plate attached to the model base. The pressures on the splitter plate show the 
recompression from the measured base pressure (an average of the nearly constant pres- 
sure indicated at the five static-pressure orifices without the splitter plate) in the wake 
neck region followed by a pressure decay farther downstream. Also shown is an exact 
calculation by the inviscid method of characteristics for a wedge followed by a flat plate 
(no rearward step). Beyond the recompression, the wake static pressure agrees with the 
calculated pressure and thus indicates that the recompression region has only a local 
effect on the inviscid-wake development. The greater displacement effect which would 
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occur On the splitter plate than in the wake must produce the higher pressures seen on 
the plate. One would expect the inviscid calculation to represent a lower bound on the 
static pressure, and certainly the splitter plate would be an upper bound. The close 
agreement of the measured and calculated wake static pressures gives further credence 
to the use of the direct measurements, as discussed earlier, for subsequent calculations 
of the viscous -wake flow properties. 

The wake flow properties were computed from the survey data as follows. First, 
the static pressure was assumed to be constant across the viscous wake. (See ref. 8.) 
From the measured static and pitot pressure the local Mach number is found from 
normal-shock relations, and the velocity and static-temperature ratios were found by 
using (ref. 9) 
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Velocity and temperature ratios for x/t = 6.0 are given in table I along with the wake- 
edge conditions. 


Transition Detection 

The occurrence of transition and the existence of turbulent flow in the wake were 
determined by surveys with a constant -cur rent hot-wire anemometer. The hot-wire rms 
voltage fluctuations across the viscous wake are shown in figure 8(a). In each profile the 
wake recompression shock appears as the region of large rms voltage outside the viscous 
wake; the locations of the high fluctuations agree with the location of the wake recompres- 
sion shock found from schlieren photographs and shown in the figure. At x/t = 1, the 
hot-wire data show little turbulence in the viscous wake, but considerable turbulence 
develops between x/t = 2 and x/t = 6. Somewhere upstream of x/t = 4, an apparent 
increase of the turbulence (above the background level) at the middle of the wake is 
observed; according to the criterion used by Demetriades (ref. 10), an increase in the 
turbulence at the middle of the wake is indicative of transition to turbulent flow. At 
x/t = 6, the rms voltage profile shows considerable turbulence throughout the viscous 
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wake. Therefore, transition seems to occur at x/t '=4. A further check on the transi- 
tion location uses a criterion put forth by Demetriades (ref. 11): that wake transition 
can be observed as a change in the x variation of the square of the reciprocal of the max- 
imum velocity defect, that is, the defect at y = 0. This criterion also shows transition 
to occur around x/t == 4. (See fig. 8(b).) The wake transition Reynolds number 

PpU„ . 

would then be about 10^ based on local flow conditions and assuming the "sticking dis- 
tance" ^the high Reynolds number upstream limit of the transition point Xg^j^tj to be 
about 1. 

Demetriades' measurements of transition locations in the wakes of wedges at a 
free-stream Mach number of 6 (ref. 10) may be compared with this transition mea- 
surement. A local Reynolds number is more appropriate to compare the results of 
Demetriades for a free-stream Mach number of 6 with the present data for a free-stream 
Mach number of 20.4, There are two reasons for this: first, Batt (ref. 8), using local 
instead of free-stream conditions at M = 6 on a 20° total-angle wedge, has found only 
a slight change in the transition Reynolds number, and second, the wake -edge Mach num- 
ber for both the Demetriades data and the present data would be about 6. From refer- 
ence 10 the transition sticking distance Xsd would be about 1.0 and the transition 
Reynolds number Rtr would be about 4 x 10^, somewhat higher than observed herein. 
The lower transition Reynolds number in the present study could be due to the "unit 
Reynolds effect" observed in reference 12 and common to wind-tunnel boundary-layer 
transition studies (ref. 9), that is, the dependence of the transition location upon the facil- 
ity free-stream disturbance level. In hypersonic tunnels the dominant source of free- 
stream disturbances is sound radiated by the turbulent nozzle -wall boundary layers 
(refs. 9 and 13). The sound level in the present tests would be much greater than in the 
lower Mach number studies owii^ to the rapid increase in sound intensity radiated from 
a turbulent boundary layer as Mach number increases (ref. 13). 

Theoretical Predictions and Comparisons With Experiment 

Theoretical calculations . - To compute the wake development from the initial input 
profiles, an implicit finite -difference method was used to solve the equations of motion. 
The method used was that developed and programed by Sinha, Fox, and Weinberger 
(ref. 14) but modified to account for pressure-gradient effects correctly. For flows hav- 
ing a streamwise pressure gradient, the continuity equation as formulated and programed 
in reference 14 is not satisfied by the nondimensional Von Mises transformation used in 
the analysis therein. (This becomes apparent when the continuity equation is nondimen- 
sionalized and expanded.) It should be pointed out that the same improper formulation is 
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used by Zakkay and Fox in reference 15. The proper formulation is given in the follow- 
ing equations. Only single-component, ideal gases are considered herein. The equa- 
tions of motion are 

Continuity: 1 


|-(puyi)*^(pvyi)=0 


Momentum: 


9u 

pu — + pv 
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9u ^Pe __-j 9 
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Energy: 


9T 
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where j = 0 or 1 for two-dimensional or axisymmetric flow, respectively. 


(3) 


By introducing a Von Mises transformation 


(1 + j)TiTy = puyi 
(1 ,+ = -pvyj 


(4) 


the continuity equation is satisfied, and the momentum and energy equations are 
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Nondimensional variables can now be introduced: 


; .y-2 1 dx 
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^b(Pe'^e) 




( 6 ) 
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and 





P 



C =-BL- 



where Xj is the initial value of x, Rg = Pg^e^b/^e’ ^ characteristic body 

dimension. The momentum and energy equations become 



where the subscripts t and f denote differentiation with respect to t and respec- 
tively. The initial and boundary conditions for the wake are 

u(0,f) = iii(T) f (0,T) = Ti(T) 

u-(Lo) = T-(?,0) = 0 lim = lim T(It) = 1 

T_»oo T— 00 


At this point, these equations and the dependent and independent variables are now 
the same as those in reference 14 except for the additional underlined terms and the fact 
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that T- is not related to ^ as indicated in reference 14. Note that for flows with 

A 

zero pressure gradient the underlined terms are identically zero, and in this case the 
analysis of reference 14 gives correct results. The modification of the program of ref- 
erence 14 for the correct computation of flows with pressure gradients is quite straight- 
forward; this modification was made by simply including the underlined terms in the trans- 
formed equations at the appropriate locations in the program. (Similar corrections for 
the species equation for flows with dissimilar gas mixing would occur.) 

Turbulence models . - Equations (7) and (8) are applied to either laminar, turbulent, 
or intermittently turbulent flow by using transport models for the transport coefficients 
p and a. The turbulent-flux terms in the mean-flow equation are modeled by an eddy 
viscosity and a turbulent Prandtl number; that is. 


Turbulent shear stress = /ij, 



Turbulent heat flux = 


^T 8h 
ffrj, dy 


} 



( 9 ) 


Where a,j, is the turbulent Prandtl number. The intermittent nature of the viscous wake 
flow is accounted for by an intermittency factor yrp so that the total shear stress (and 
correspondingly for the heat transfer) is 


Total shear stress = u ^ = 

9y 



(10) 


where yj, is 1 for fully turbulent flow and 0 for fully laminar flow. Since a = Cp/j/k, 
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T Mi 


1 + yrp 


Mfp 

Mi 


For simplicity, the laminar Prandtl number was assumed to be equal to the turbulent 
Prandtl number a,p, so that o = Oj,. An approximation made by Townsend and given by 
Hinze in reference 16 for the distribution of the intermittency was assumed herein for all 
calculations: 


yt - 



( 11 ) 
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where yj/2 velocity -defect half -width of the wake ^.e., where u(y) - uq = 

The constant Cy was taken as 0.142, a value which gives close agreement with the inter- 
mittency distribution given in reference 17. Demetriades (ref. 18) has measured the vari- 
ation of y-p in the wake of a thin plate at a Mach number of 3 and has found the same 
dependence of y-p upon y as in equation (11) and nearly the same value of Cy. 

Different formulations for the eddy viscosity (representative of those found in the 
literature) were examined. Calculations were made by assuming an eddy viscosity pro- 
portional to the velocity defect in the wake; that is, 



Mt ~ ^uPeyi/2^^6 ~ 


( 12 ) 


or, proportional to an integral of the mass-flow defect, 


Mt “ ^m 



(PgUe - pu)dy 


(13) 


The first formulation is the classical Prandtl hypothesis (ref. 19) and the latter has been 
suggested by Schetz (ref. 20). Both methods predict a constant eddy viscosity across the 
wake. 

In addition to these models, calculations were made on the basis of a mixing-length 
type of formulation for the eddy viscosity 


p.>p = pZ^ 




(14) 


which has been suggested in the literature by Prandtl (ref. 19). The yelocity-gradient 
term and the mixing length Z are the usual quantities introduced in the typical approach 
for turbulent boundary layers, but the term proportional to the second derivative of the 
velocity was suggested by Prandtl to modify the usual mixing-length model which predicts 
vanishing eddy viscosity at the wake plane of symmetry, a result inconsistent with obser- 
vations. In all calculations herein, it was assumed that ^zyyjy-2j = 0.32 and Z = Z^. 

Since yj^/2 '"I ^ wake width), the variation of Z chosen is nearly Z -0.086, which 
is a good approximation for mixing lengths in the outer part of turbulent boundary layers 
(ref. 21). This value of is close to the value obtained by Schlichting in low-speed 

wakes (ref. 19). 

Comparisons of predictions and experimental data .- Initial velocity and temperature 
profiles for the finite-difference calculations were obtained from the experimental data at 
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x/t = 6.0. The downstream development of the wake was computed using each of the tur- 
bulent transport models discussed previously, and one calculation was made assuming 
laminar transport. For the laminar transport properties, the assumptions a = 1.0 and 
Pju = PqIJLq were made. Experimental and predicted velocity and temperature defects at 
the plane of symmetry (the maximum defects) are shown in figures 9(a) and 9(b). Edge 
conditions were obtained from the values computed at the edge of the viscous wake or just 
downstream of the wake shock for those stations where the viscous wake edge is not well 
defined (x/t less than about 5). The departure of the experimental data from the laminar 
calculation beyond the initial station is further evidence that transition to turbulent flow 
occurs in the wake upstream of x/t = 6.0. 


Two calculations shown in figures 9(a) and 9(b) were performed with the modified 
mixing-length model: one calculation with a,p = a^=1.0 and the other with o,p = a^ = 0.68. 
The latter value equals the actual molecular Prandtl number for helium at the tempera- 
tures present in the wake. The Prandtl number assumption appears to have a measurable 
effect upon the temperature-defect variation but only a slight effect on the velocity defect. 
Although better results would probably occur for a value of Orj, between 0.68 and 1.0, 
either of the present calculations is considered to yield good predictions. 


Calculations were made with the classical mixing-length theory (Zj = o) and assum- 
ing still that ^i/yi/ 2 ) = 0.32, but the results were unrealistic. (See figs. 9(a) and 9(b).) 
Although the turbulent viscosity was as much as 30 times greater than the laminar vis- 
cosity, the calculations gave almost as slow relaxation (i.e., less mixing) of the velocity 
and temperature, defects as that occurring in the laminar solution. This result was not 
expected, since the primary effect of the additional term in the modified mixing-length 
theory is to give nonzero values for the eddy viscosity at the plane of symmetry; else- 
where the eddy viscosity is not greatly changed. For the mixing-length theory, the rapid 
decrease of the eddy viscosity as y approaches zero apparently is difficult to accom- 
modate in the solution. Schlichting (ref. 19) has pointed out that analyses of low-speed 
free shear layers using the classical mixing-length theory have the general property that 
discontinuities in Uyy occur in the solutions; for example Uyy becomes infinitely large 
at y = 0 in a wake. The unrealistic results obtained herein may be a manifestation of 
such a singularity. Schlichting (ref. 19) also states that the modified mixing-length theory 
is free of such difficulties. 


Through appropriate choices of the constants in the eddy-viscosity formulations of 
the algebraically simpler models, Prandtl's hypothesis (velocity-defect model) and Schetz' 
model, good predictions of the maximum velocity and temperature defects could be 
obtained. (See figs. 9(a) and 9(b).) In the velocity-defect model, the value of Ky used 
was 0.036. Lees and Hromas proposed a velocity-defect model with Ku = 0.040 and a 
transformed wake thickness (ref. 22) instead of the velocity-defect half -width y\/2 'ised 
herein. For the density ratios present in these experiments, these changes in the model 
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would not have large effects on the eddy -viscosity levels. In Schetz' model the value of 
Kj^ was 0.005; this value is an order of magnitude lower than the value recommended by 
Schetz, Kjjj = 0.060. (In an evaluation of axisymmetric jet mixing, Eggers (ref. 23) has 
similarly found Schetz' recommended constants too large by a factor of 6.) This model 
gave nearly the same solution as the velocity-defect model and for this reason will be 
omitted in later comparisons of theory and data. 

Velocity and temperature distributions from the experiments and the computations 
are shown in figures 10(a) and 10(b). All the theoretical models gave solutions which 
yield good predictions of the velocity distributions. Because of the nearly identical veloc- 
ity distributions, only the modified mixing-length theory (a-j. = l.O) is shown in figure 10(a). 
The temperature distributions (fig. 10(b)) do show some differences between the calcula- 
tions. Again, as for the maximum temperature defects, the two calculations by the modi- 
fied mixing-length theory indicate that improved calculations of the temperature distribu- 
tions near y = 0 can be obtained for. Orj> between the values shown. The velocity -defect 
model, although yielding good predictions of the maximum temperature defects, predicts 
temperature profiles which are slightly fuller than the experimental profiles as do the 
modified -mixing -length calculations. Nevertheless, all the calculations are considered to 
give good predictions. 

The momentum thickness 9 and displacement thickness 6* of the wake are shown 
in figure 11. Also shown are the velocity-defect half-widths The increase in 

momentum thickness with x/t occurs because of the pressure gradient in the wake. To 
check the consistency of the data, the initial value of 6 and the experimental values of 
6*, Me, and dpg/dx were used in the integral momentum equation (ref. 19) to calculate 
the momentum-thickness growth. The calculated variation (see fig. 11) agrees well with 
the experimental variation of 6 with x. The generally good agreement of the data and 
the results of the finite -difference solutions would be expected in view of the close profile 
predictions previously discussed. 

Turbulent transport models which give the eddy viscosity proportional to a velocity 
or mass-flow defect, such as the Prandtl hypothesis or Schetz' model, would be expected 
to be appropriate for a wake with near self-preserving similarity. Townsend (see ref. 18) 
has shown that at low speed a self -preserving similarity can be attained in zero-pressure- 
gradient wakes for Reynolds numbers, based on the wake momentum thickness, greater 
than 400. Demetriades has observed in a two-dimensional supersonic wake (ref. 18) that 
the self -preserving similarity of incompressible far wakes is also reached in compressible 
wakes. The Reynolds number based on momentum thickness at the initial station herein 
is Rq ~ 800. The question exists whether the magnitude of this Reynolds number is such 
that a self-preserving similarity can be established in the present wake with nonzero pres- 
sure gradient. It was of interest to determine whether similarity was approached in the 
present experiments. 


15 



An analysis of the admissible similar solutions, or the appropriate similarity vari- 
ables, for compressible wakes with a pressure gradient (as described by eqs. (7) and (8)) 
was not attempted. Instead, similarity was investigated by using the similarity variables 
firmly established in the literature for wakes with zero pressure gradient (refs. 16 to 18) 
in the belief that the streamwise pressure gradient would have only a small effect on the 
appropriate variables. In figure 12 this approach is taken, and the velocity and tempera- 
ture defects in the wake are shown divided by the corresponding maximum defects and 
plotted against the transverse coordinate y scaled by the half-width of the wake velocity 
defect y-1^2' velocity-defect data display "similarity" and the correlation agrees 
closely with the results obtained by Townsend at low speeds. The temperature-defect 
data also display similarity even though the transverse coordinate is scaled by the 
velocity-defect half-width. (Note that the low-speed velocity defect is also shown in the 
temperature -defect correlation for comparison.) The degree of similarity shown is 
rather surprising in view of the results obtained by Demetriades (ref. 18), who finds that 
at Me = 3 the wake of a flat plate does not attain self-preserving similarity until the 
distance from the virtual origin of the turbulence (see fig. 8) exceeds 10000. Herein the 
minimum value of 0 is about 0.02t, and therefore, similarity would not be expected 
upstream of x/t ~ 20. 

Although the present experiments contain pressure-gradient effects, it is of interest 
to determine the turbulent Reynolds number (ref. 18) defined for wake flows with zero 
pressure gradient as 



where x^q is the virtual origin of the turbulent wake. Using the linear curve fit down- 
stream of transition in figure 8 and an average value of 0 gives R-p = 12.0, which is 
comparable to the Demetriades (ref. 18) value, Rx = 13.0. 

C ONC LUDING RE MARKS 

Experiments have been performed to determine the flow in the turbulent wake pro- 
duced by a two-dimensional wedge in hypersonic helium flow. From these experiments, 
mean-flow profiles were obtained and compared with finite-difference calculations of the 
subsequent wake development with various models inserted for the turbulent transport. 
The comparison of calculations and experimental data shows that by proper choice of 
empirical constants, the turbulent transport models that employ eddy -viscosity formula- 
tions based on velocity or mass-flow defects can yield accurate predictions of hypersonic, 
turbulent wake flow including the effects of pressure gradients. However, the classical 
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mixLng-length theory was found to yield unrealistic results when used in the finite - 
difference calculations. When the classical mixing-length theory is modified to include 
a second-derivative term in the eddy -viscosity formulation, the difficulties of the classi- 
cal mixing -length theory are not encountered. In fact, the modified mixing-length theory 
yields excellent predictions of the wake development when the mixing-length constants in 
the theory are chosen to be consistent with turbulent-boundary -layer results. This link 
to the results of turbulent-boundary-layer theory is a strong point in favor of the modi- 
fied mixing-length theory over the defect-type models for which the choice of empirical 
constants is apparently flow dependent. 

The experimental data displayed similarity when velocity and temperature defects 
were scaled by the maximum defects and the transverse coordinate was scaled by the 
velocity-defect half-width. On the basis of the results of previous studies, this similar- 
ity would not be expected in the region investigated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 6, 1972. 
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TABLE I.- EXPERIMENTAL VELOCITY AND TEMPERATURE RATIOS 


x/t = 6.0; Me =5.15; 

Tt e = 93.4° C; Ue = 1850 ^; — = 9.19 


x/t = 
Tt^e = 94.7° 

= 8.98; Me=5.69; 

C; Ue = 1868S; £l = 8.57 
’ ^ s’ n 

IT OO 


x/t=10.0; Me = 5.79; 

Tt e = 93.4°C; Ue = 1872 S; — = 8.30 

Poo 

y/t 

u/ue 

T/Te 


y/t 

u/ue 

T/Te 


y/t 

u/ue 


0 

0.8524 

2.906 


0.004 

0.8968 

2.535 



0.9061 

2.448 

.011 

.8524 

2.906 


.036 

.8968 

2.535 



.9088 

2.380 

.042 

.8598 

2.816 


.071 

.9025 

2.419 



.9121 

2.317 

.073 

.8692 

2.673 


.106 

.9076 

2.319 



.9176 

2.186 

.105 

.8827, 



.140 . 

.9168 

2.177 


.180 

.9325 


.136 

.8975 

2.287 


.180 

.9274 

2.028 


.230 

.9439 

mgm 

.163 

.9138 

2.048 


.240 

.9516 

1.679 



.9601 

119 

.192 

.9269 

1.867 


.320 

.9744 

1.404 



.9799 

WSSM 

.225 

.9445 

1.678 


.390 

.9875 

1.211 



.9898 

1.195 

.253 

.9614 

1.458 


.480 

.9947 




.9951 

1.109 

.320 

.9820 

1.252 


.560 

.9968 




.9973 


.380 

.9948 

1.090 


.670 

.9981 



.740 

.9990 


.410 

.9966 

1.059 


.790 

.9987 

1.009 



.9993 


.450 

.9984 



.900 

1.0000 

1.000 



1.0000 


.510 

1.0000 

1.000 







1.0000 



x/t = 12.0; Me = 5.95; 

Tt e = 91.2° C; Ue=1874— ; ^ = 7.69 
> e s ’ p„ 


x/t=14.0; Me = 6. 11; 

Tt,e = 93.40 C; Ue=1877 f; §^=7.03 

^ oo 


x/t= 16.94; Me = 6.59; 

Tt 6 = 93.2° C; Ue = 1890S; ^=6.00 
’ ® Poo 

y/t 

u/ue 

— 

T/Te 


y/t 

u/ue 

T/Te 


y/t 

u/ue 

T/Te 

0 

0.9177 

2.295 


0 

0.9219 

2.302 


0 

0.9300 

2.120 

.010 

.9177 

2.290 


.010 

.9250 

2.300 


.013 

.9379 

2.116 

.030 

.9180 

2.280 


.040 

.9230 

2.240 


.038 

.9390 

2.100 

.084 

.9208 

2.225 


.105 

.9272 

2.204 


.106 

.9404 

2.050 

.170 

.9290 

2.030 


.196 

.9348 

2.044 


.201 

.9447 

1.927 

.260 
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1.764 


.260 

.9427 

1.868 


.280 

.9516 

1.788 

.350 

.9649 

1.493 


.380 

.9624 

1.548 


.360 

.9609 

1.604 

.440 

.9821 

1.274 


.490 

.9813 

1.306 


.450 

.9724 

1.441 

.540 

.9934 

1.119 


.640 

.9943 

1.122 


.540 

.9825 

1.297 

.640 

.9972 

1.066 


.810 

.9954 

1.060 


.630 

.9901 

1.191 

.740 

.9981 

1.049 


.980 

.9982 

1.038 


.730 

.9950 

1.128 

.840 

.9988 

1.029 


1.130 

.9990 

1.015 


.830 

.9970 

1.082 

.940 

.9998 

1.010 


1.330 

1.0000 

1.009 


.930 

.9980 

1.060 

1.040 

1.0000 

1.000 


1.450 
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1.000 


1.030 

.9988 

1.035 









1.140 

.9993 
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1.240 

.9996 

1.009 









1.360 

1.0000 

1.000 
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Figure 5.- Total-temperature distributions in wake of 22.5° half-angle wedge. Moo = 20.4; R^o = 3.4 x 10 
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Figure 6.- Lines of constant pitot pressure. Moo = 20.4; Roo = 3.4 x 10‘ 
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Figure 7.- Static -pressure distribution in wake and on splitter plate. 
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(b) Velocity defects. 

Figure 8.- Transition detection in wake of 22.5° half-angle wedge. 
Moo = 20.4; R«, = 3.4xl0^ 
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(b) Maximum temperature defects. 
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Data 
a= 1.0 
cr = 0.68 


Modified mixing-length theory 



(b) Temperature profiles. 
Figure 10.- Concluded. 
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Figure 11.- Wake integral properties. 
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Figure 12 .- Correlatious of velocity and temperature defects 
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